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! Abstract: Taking the isotropic limit A ^ 1 in a recent representation theoretic con- 

^sO \ struction of Baxter's Q-operators for the XXZ model with quasi-periodic boundary con- 
ditions we obtain new results for the XXX model. We show that quasi-periodic boundary 

^ ' conditions are needed to ensure convergence of the Q-operator construction and derive a 

■ quantum Wronskian relation which implies two different sets of Bcthc ansatz equations, 

. one above the other below the "equator" of total spin S^' ~ 0. We discuss the limit 
to periodic boundary conditions at the end and explain how this construction might be 

lO ! useful in the context of correlation functions on the infinite lattice. We also identify a 

' special subclass of solutions to the quantum Wronskian for chains up to a length of 10 



sites and possibly higher. 



1. Introduction 

Historically Baxter's Q-operator was introduced as substitute method for the coordinate 
Bethe ansatz in solving the eight-vertex model [l|, but has more recently seen wider 

applications in the field of integrable systems e.g. |@, ^, ^ 0) HHl making it an important 
and almost universal tool. To motivate his technique Baxter first discussed the concept of 
the Q-operator in the context of the six-vertex or XXZ model, where a direct comparison 
with the coordinate Bethe ansatz can be made. While our primary interest in this article 
will be the XXX model it is helpful to consider first the anisotropic or XXZ case. Denote 
by t the transfer matrix then the Q-operator is implicitly defined through the functional 
equation 

t{u)Q{u) = x{u - \Y'Q{u + 1) + x{u + D'^Qiu - 1) (1.1) 

where x is an explicitly known function (the quantum determinant) and M the number 
of lattice columns, respectively the number of sites in the spin-chain. In addition, to 
this relation, known as TQ equation, one usually requires a number of properties such 
as "analyticity" of the Q-operator in the spectral variable u and that [T(n),Q(n')] = 
[Q{u)^Q{u')] = for an arbitrary pair u,u' G C. The latter commutation relations allow 
one to discuss the TQ equation on the level of eigenvalues and this is where one makes 
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contact with the coordinate Bethe ansatz |10| which determines the spectrum of the transfer 



matrix in terms of the solutions {vk}^^i to the Bethe ansatz equations |11, 

( "tt'ivo! )"-!! .■=i.2,...,»=M/2-s'. (1.2) 

\smh7(fj + 1/2) J fj-, smh7(fj - Vfc + ^) 

Here 7 is the crossing or couphng parameter of the six-vertex model and S"^ > the total 
spin-operator. Postulating that the eigenvalues of the Q-operator are of the form Q 

j=i smh 7 

the outcome of the coordinate Bethe ansatz then implies the TQ relation ( |1.1| ), which is 
the starting point for the construction of the operator Q. Note that this line of argument is 
based on the essential assumption that the coordinate Bethe ansatz yields a complete set of 
eigenstates of the transfer matrix with a finite set of Bethe roots Vj. It is this assumption, 
which has to be treated with care in the isotropic limit 7—^0 yielding the XXX model^. 

The transfer matrix of the XXX model as well as the associated Heisenberg spin-chain 
are SI2 symmetric, whence their eigenspaces decompose into SI2 modules. As is well known 



the finite solutions to the XXX Bethe ansatz equations (first derived by Bethe in |10| albeit 
in a different form), 

Vj-i/l V' _ v--v^-i ^ ^^^^^ 

now only yield the highest weight vectors in each s/2 module The remaining states 
within each module are obtained through the action of the symmetry algebra and have 
been referred to as "non-regular" Bethe states as they involve "infinite rapidities" in the 
particular parametrization used in (|1.4|); see e.g. |17| for a discussion how to recover the 



non-regular Bethe states through a limiting procedure. Thus, the obvious ansatz 

n 

Q(.u)=l\(.u-Vj) (1.5) 

for the eigenvalues of an XXX Q-operator becomes problematic due to the presence of 
"infinite rapidities" , or more precisely not all states correspond to finite solutions of the 
Bethe ansatz equations ( |1.4| ) . Clearly, there are ways out of this dilemma, either by choosing 
a different parametrization such that all rapidities stay finite (this is for instance the case 
in the coordinate Bethe ansatz, where the non-regular Bethe states correspond to the case 
that multiple quasi-momenta vanish), or by (continuously) breaking the SI2 symmetry in 
such a manner that the assumption on the completeness of the Bethe ansatz becomes 
applicable again. 

In this work we shall do the latter by introducing quasi-periodic boundary conditions. 



see e.g. ||l^, |T^, |2^, |21[]. This has the advantage that all relevant algebraic properties 



^Similar problems occur for the XXZ model at roots of unity |l^ due to a partial loop algebra 
symmetry [|l5||. 
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needed for the quantum inverse scattering method [22| stay intact and that we can take 
at the very end the hmit to periodic boundary conditions making contact with previous 
investigations of Q-operators for the XXX model. Of particular interest will be aspects 
which are not accessible through the coordinate Bethe ansatz, namely the existence of two 
linearly independent solutions, say Q^, to the TQ equation and, closely related with this 
question, the derivation of the following quantum Wronskian identity 

ioQiiu - ^)Q-{u + i) - ^-iQ+(n + §)g- (n - 



to — UJ 



Xiu) , (1.6) 



which is a new result. Here uj = exp{i(j)) is the twist parameter associated with the quasi- 
periodic boundary conditions and x is the aforementioned quantum determinant, but now 
of the XXX model. Since the latter is explicitly known, e.g. x(^) = for the homoge- 
neous case, one can employ the quantum Wronskian ( |l.6D rather than the generalization of 
the Bethe ansatz equations ( |l.4D to twisted boundary conditions when solving the model. 
Namely, making the ansatz (which will be justified through our construction of in the 
text) 

n M—n \J 

QUu) =\{{n- vj) and Q^(n) = H (^ " vj\ n = - - (1.7) 
i=i i=i ^ 

for the eigenvalues of the two solutions to the TQ equation, the roots = v^{uj) are 
determined through ( |1.6| ). Here denotes the total spin component in the direction singled 
out by the quasi-periodic boundary conditions. Note that upon setting u = v^+i/2, —i/2 



the identity (1.6) implies two different sets of Bethe ansatz equations, one above, the other 



one below the equator = 0. Due to the quasi-periodic boundary conditions, w 7^ 1, the 
Bethe roots are all finite and the number of solutions matches the dimension of each 



fixed spin-sector signaling completeness; compare for example with the discussion in ||25 

As discussed above this ceases to be true in the limit a; — > 1 corresponding to periodic 
boundary conditions. From (|l]^) we infer that this limit might indeed be singular unless 
the numerator and denominator vanish simultaneously. We will compare the outcome of 
this article with the findings for periodic boundary conditions by Pronko and Stroganov 



1 24], who have presented a similar quantum Wronskian without the denominator at a; = 1 
and a different degree for the second solution Q~ , namely degQ" = M — n + 1. Their 
Wronskian relation can be numerically solved but the resulting number of solutions is in 
general much smaller then the dimension of the state space (™). In light of the previous 
remarks on the s/2 symmetry this is not surprising as their solutions only yield the highest 
weight vectors in each module. Taking the limit w ^ 1 in the explicit solutions to ( |l.(i| ) 
for small chains we indeed find that of those solutions which stay finite, both approach 
the solution of Pronko and Stroganov. We shall comment on this in more detail in the 
text, see section 5.2. 

The appearance of singularities in the limit of periodic boundary conditions can also 
be understood from the explicit construction of the Q-operator for twisted boundary con- 
ditions. The latter is given as the trace of a monodromy matrix with infinite-dimensional 
auxiliary space. In order to obtain a well-defined object one must ensure convergence 
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of the trace. As we will see in the text this actually requires the introduction of quasi- 
periodic boundary conditions. Previous constructions of Q-operators for the XXX spin- 



chain 27 1 have been for periodic boundary conditions only, where Q has been 

represented as an integral kernel (see also for a related XXZ construction). 

In contrast the limit of the transfer matrix from quasi-periodic to periodic boundary 
conditions is well-defined. In fact, this applies to all higher spin transfer matrices which 
can be expressed in terms of as follows. 



(u - ^] 



tin; x) = hm ^^^^(^ " f + f ) " ^"^Q^^^^ + f ^' . 

-I . . . . 1 



+ 1 — LO 



When X = n €z N>o the function t{u; x = n) gives the spectrum of the transfer matrix 
with spin s = (n — l)/2 in the auxiliary space. However, if we take x to be an arbitrary 
complex parameter, the resulting spectrum belongs to a generalized transfer matrix used 



in the discussion of correlation functions for the infinite chain |29, 30, 31, 32 1. This result 
is the analogue of a previous discussion for the XXZ model [33, 34] and the discussion 
presented here is in accordance with these earlier results for the more general case when 
7 7^ 0. At the moment there appears to be no construction of an Q-operator for lo = 1 
which allows to define ( |1.8| ) . This is one of the main reasons for the construction presented 
in this paper. 

In section 2 the basic definitions of the XXX model and its fusion hierarchy is stated. 
Section 3 contains the construction of the Q-operator which is simply the isotropic limit 



(7 — > 0) of earlier constructions for the XXZ model [34|. We briefly address the afore- 
mentioned conditions for convergence due to an infinite-dimensional auxiliary space and 
state the relevant functional equations with the transfer matrix. We omit most proofs for 
those results which readily follow from taking the isotropic limit in the XXZ construction. 
For instance, the eigenvalues of the Q-operator are discussed by making contact with the 
algebraic Bethe ansatz discussion in |^3|. By comparison with the analogous results for 
the XXZ model it is shown that the Q-operator factorizes into two linearly independent 
solutions to Baxter's TQ-equation. We discuss how they are related via spin-reversal. The 
relation with the fusion hierarchy and its analytic continuation ( |1.8| ) to "complex spin" is 
presented in section 4. Section 5 gives the quantum Wronskian relation between the two 
independent solutions to Baxter's TQ equation, which is then compared against the one of 



Pronko and Stroganov [24|. A special subset of solutions to the twisted quantum Wronskian 
( |1.6| ) is also discussed based on numerical results for chains of even length < 10. Their 
associated Bethe roots obey identities which imply (and are therefore more fundamental 
than) the Bethe ansatz equations. The conclusions are stated in section 6. 



2. Definitions 

Let us start by introducing our conventions for the definition of the XXX model. Denote by 
{(j^ = cr^, cj^ = fj^, = cr^} the Pauli matrices acting on and let P be the permutation 
operator, P(z; ®w) = w ®v. Then the basic ingredient for constructing the XXX model is 
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the following simple solution to the Yang-Baxter equation 

r(A) = A + i+ ct"® a" = A + Pg End(C2 ® C^) (2.1) 

a=l 

Note that here we have changed our conventions from that in the introduction as it sim- 
plifies some of the following computations. Another definition of the XXX r-matrix is also 
commonly used in the literature, 

3 

f{u) := ir{-iu - 1/2) = u + i ^ a'^ (E) a'^ = u - i/2 + iF . (2.2) 

a=l 

Both definitions only differ by a re-parametrization of the spectral parameter, A —iu — 
1/2, and an overall factor i = The functional relations and equations stated in the 



introduction refer to this last convention (2.2) 



In terms of (^]^) the transfer matrix of the inhomogeneous XXX model with quasi- 
periodic boundary conditions is defined as follows, 

t^(A) = TrL^-^®VM(A-AA/)---ri(A-Ai) , cj = e'^ . (2.3) 

Here the trace is taken in the first factor of the r-matrix, i.e. G End(C^)®^^. The 
set {Xm}m=i some arbitrary generic inhomogeneity parameters, while the parameter 
CO = exp(i(/)) incorporates the twist angle (p which for the moment is allowed to be a 
generic complex number, but can be specialized later on to real values in order to ensure 
hermiticity. In the homogeneous limit Ai = ... = Aj\/ = its meaning becomes apparent 
when writing down the associated spin-chain Hamiltonian 



dX (A + 1)^^ 
with the boundary conditions 



A=o 2 



I i^m ■ ^m+i - 1) (2.4) 



a^+i ± ialj^^ = uj^^af ± ia^) and = af . (2.5) 

These boundary conditions break for uj ^ 1 the spherical symmetry of the Hamiltonian 
which unlike in the case of periodic boundary conditions is not s/2 invariant. However, 
there is an axial symmetry, i.e. the total spin operator 

1 ^ 

5' = 2 E < (2-6) 

m=l 

is preserved. This breaking of the spherical symmetry is significant for the Bethe ansatz 
analysis of the spectrum as for quasi-periodic boundary conditions all eigenvectors become 
regular Bethe states. In the case of periodic boundary conditions this is only true for the 
highest weight state in each s/2-iiiodule spanning one of the degenerate subspaces of the 
transfer matrix respectively the Hamiltonian. This fact also plays an important role in the 
construction of the Q-operator. 
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Besides the transfer matrix and the Hamihonian it wih be convenient to discuss the 
entire fusion hierarchy of the XXX model. To this end consider the ChevaUey-Serre 
generators of sh, 

[^,e]=2e, [/!,/]= -2/ and [e,f] = h, (2.7) 

then the fohowing defines a well-known Verma module i^x depending on a complex param- 
eter X G C, 

7r^(e) \k) = {x- k)k \k-l), 7ra.(e) |0) = (2.8) 
'^x{f)\k) = + 

'Kx{h)\k) = {x-2k-l)\k) , A; = 0,l,...,oo . 

It is this Verma module which will form the auxiliary space for the Q-operator. Note that 
if X = ri G N>o and one invokes the truncation condition iTxif) \n) = 0, the ra-dimensional 
subspace spanned by the vectors {|fe)}^~Q gives rise to the finite-dimensional modules 
ttC"-!) known as spin s = (n — l)/2 representations in the physics literature. Set 

/ A ^ f \ 
^^^'^ " ( e ^ A - ^ ) ^ ^ ^^'^'^^ ^^'^^ 

then 

ii2(A)Li3(A + A0r23(A0 = r23(A')Li3(A + A')Li2(A) (2.10) 

(n) 

and the higher spin transfer matrix tl, is defined through 

4")(A)= Tru^^'LM{X-XM)---Li{X-Xi) . (2.11) 

7r(") 

The two distinguished elements in this hierarchy are the previously introduced transfer 
matrix t^^ = t^^ of spin 1/2 and the quantum determinant % corresponding to the trivial 
representation of spin 0, 

x(A)=t(0)(A)= n(A-A„ + i). (2.12) 

m=l 

From these two elements all the other members of the fusion hierarchy can be generated 
via the the functional equation 

4") (A + ^^)tj^\X) = t(°)(A + ^)4""^')(A + i) +i^°HA - ^)4"-')(A + ^) . (2.13) 

Instead of solving this functional relation in terms of t!^\ t^^\ which leads to quite involved 
formulae, it is simpler to consider an auxiliary linear problem, Baxter's TQ equation, which 
we discuss next. 
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3. The Q-operator and its spectrum 

We extend the definition of tlie higher spin transfer matrix to the infinite-dimensional 
Verma module (^) introduced above and set 

g^(A;x) = Tra;'^®%/(A-AM + f)---Li(A-Ai + f) . (3.1) 

This definition of the Q-operator coincides with the isotropic limit of the definition for the 
XXZ model |34]. Note that the trace runs now over an (half) infinite-dimensional space, 
whence it is crucial to have quasi-periodic boundary conditions which upon the right choice 
of the twist parameter u; ensure convergence [33|. 

Since the matrix Qij(A;x) preserves the total spin, [Q^{X; x), S^] = 0, its matrix ele- 
ments do always contain the same number of the Chevalley-Serre generators e and /. Using 
the Casimir relation, 

^x(C) = ^^, C = h^/2 + h + 2fe, (3.2) 
we deduce that it suffices to ensure that the following expressions are finite 

oo 

TT{uj^^^h"'}=uj'' Z^'^^~\x-2k-ir <oo, m = 0,1,2,..., M . (3.3) 

"^^ k=0 

This is obviously guaranteed as long as > 1. Employing the geometric series to compute 
the trace, we then analytically continue this operator from the region of convergence to the 
whole complex cj-plane. Note that there remains a pole at w = 1. 

For instance, by construction Qcj(A;x) is a polynomial of degree M in A and we have 
for the coefficient of the highest power A*^, 

oo 3, 

QUX;x) = TV{u;n A^^ + ... = y uj^-^'^-^ + ... = A*^ + ... , (3.4) 

k=0 

where the last expression can be continued with respect to uj from the region of conver- 
gence into the complex plane. Henceforth, this analytic continuation from the region of 
convergence shall always be implicitly understood. 

The crucial property of the Q-operator is the following functional equation, which 
strictly speaking is not yet Baxter's TQ equation, 

M M 

t^(A)Q^(A;x) = Q^(A + 1;2;-1) J] (A - A™) + Q^(A - 1; 2; + 1) n (A - A„ + 1) . (3.5) 

m=l m=l 

We omit the proof as it follows from taking the isotropic limit in the analogous XXZ 
relations; see e.g. 34]. The difference with Baxter's TQ equation is the fact that the 

additional complex parameter x originating from the definition of the Verma module also 
shifts, instead of only a shift in the spectral variable A. Thus, the above equation should 
rather be seen as an extension of the fusion hierarchy to "infinite" spin. Nevertheless, 



-7- 



A Q-operator for the twisted XXX model 



the solutions to Baxter's TQ equation are obtained from Q(^(A;x) through special limits. 
Namely, as we will discuss below we have the following factorization 

Q^{X;x) = ^^Q+{X)QZ{X + x), (3.6) 

U! — UJ ^ 

where are two linearly independent solutions to Baxter's TQ equation 

M M 

t.WQtW = ^^^Qti^ + 1) n (A - A^) + iv^^QtiX - 1) n (A - A^ + 1) . (3.7) 

m=l m=l 

We now turn to the discussion of the spectrum of the Q operator where we will explain in 
more detail the above factorization into the solutions Q^. 

3.1 The algebraic Bethe ansatz analysis of Q 

In the context of the XXZ model the spectrum of the Q-operator has been analyzed p3|] 



using the formalism of the algebraic Bethe ansatz [22|. We recall that for twisted boundary 
conditions there is no problem with the Bethe ansatz as the sl2 symmetry of the XXX 
model is broken and all eigenstates of the XXX transfer matrix are regular Bethe states; 



see for instance the discussion in |23| where the completeness of the Bethe ansatz in a 
neighbourhood of w = has been shown. Using the analogous algebraic relations as in the 
XXZ case (Esl, one can show that the Bethe states are eigenvectors of the Q-operator^. 



Namely, decomposing the monodromy matrix of the XXX model in the usual manner 

(3.8) 

one considers for n = M/2 — 5"^ > an "admissible" solution p^] to the Bethe ansatz 



equations above the equator (note that according to the conventions ( |2.1| ) and (2.2) the 
corresponding Bethe roots are related by — > —iv'j — 1/2) 

n M n M 

n(c^+-c; + i) n (^+-A™)+^n(e+-e+-i) n (e+-A™ + i) = o. (3.9a) 

j=l m=l j=l m=l 

Then it follows from the Yang-Baxter equation that the matrix elements {Qki}k,i(^n of the 
monodromy matrix 

Q(A) = oj^'^^LMiX - Aa^ + f ) • • • Li(A - Ai + f ) 

with respect to the infinite-dimensional auxiliary space corresponding to Tlx satisfy certain 
commutation relations with the Yang-Baxter algebra {^4, C, D}, for example [^] 

akOii+i 

+ ^ Qu,i+iWA{0 - ^ Qk+i,i{X)D{0 + ^^Qfc+i,,+i(A)C(0, 



^At the moment this has only been carried out for Bethe states with n < 4 due to the complicated and 
numerous unwanted terms, see the appendix in p^■ However, in the case of the XXZ model alternative 
proofs (based on functional relations) exist [p6| which match the algebraic Bethe ansatz result for 
arbitrary n. The spectrum for the XXX model presented here is the isotropic limit of the XXZ result 
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where the coefficients are determined through the matrix elements of the L-operator (p.S|), 
Ofc = A — ^ + X — /c, (3f^ = l, = {x — k)k, 5k = X — ^ + k + \. (3.10) 



For a complete list of the algebraic identities we refer the reader to ^3[. Employing these 
commutation relations one can identify the eigenvalue of the Q-operator on a Bethe state. 
Denoting by |0) the pseudo vacuum, i.e. the state with all spins up, the Bethe vector 
associated with an admissible solution to (3.9a) is an eigenstate of the Q-operator with 
eigenvalue 



Q^(A; x)B{i+) . . . B{Ct) |0) = ^^Q+(A)g-(A + x)S(C+) • • • Bi^) |0) (3.11) 

UJ — Lu 



where 



«(A)= n(A-e;) (3.12) 



and 

<3:(A)^(.-.-)<3i(A)|g^^a|^^. (3.13) 

Notice that is polynomial in A due to the the Bethe ansatz equations. In fact, by the 
very construction of the Q-operator it must be a polynomial of degree M — n, 

M~n 

QZW= n (A-C7) . (3.14) 

Exploiting the completeness of the Bethe ansatz for generic quasi-periodic boundary con- 
ditions and inhomogeneity parameters [^], we obtain the factorization of the Q-operator 
into the previously introduced, linearly independent solutions Q^ of Baxter's TQ equation 
(|3.7| ). We might define them implicitly as operators through the following limits 

hm Q^(A; x) = -^^^Q+{X)Q- (0) (3.15) 

and ^ 

hm QUX;x) = ^^Q+(0)Qj(x) . (3.16) 

We shall denote the operators and eigenvalues by the same symbol. In contrast to the 
XXZ case [U] the operators Q^(0) are not easily determined and we are missing at the 
moment concrete operator expressions for them. However, explicit computation of the Q- 
operators in the various spin-sectors for small lattice sizes (M < 6) shows that the following 
expressions drastically simplify 

uj'^QU^; -A)Q.(0; 0)-i = Q+(A)Q+(0)-i (3.17) 

and 

uj-^QUO;0)-'QUO;X) = Q-(o)-iQ-(A) . (3.I8) 

Both ( 3.12 ) and ( ^.14 ) are obviously solutions to Baxter's TQ equation ( |3.7D and are 
normalized to the identity matrix at the origin A = 0. The inverse matrices exist as long as 
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none of the Bethe roots vanishes, which is the case as long as w 7^ 1. Despite this lack of 
information on the normalization constants, our Q-operator analysis yields computational 
advantages. Before we address the latter let us first discuss the relationship between Q^{X) 
under spin reversal. 

3.2 Spin reversal 

Define the spin reversal operator by setting 9^ = nm=i '^m- I-^^e to the twisted boundary 
conditions spin reversal symmetry is broken and we have for the transfer matrix the identity 

mtux)^ = t^-iW ■ (3.19) 

Let us now investigate the transformation of the Q-operator under spin reversal. From the 
equality 

/ _A - 1 -I- ^^-ti -p \ 

{l^a^)L{X){l^a^) = -i _J 2 (3.20) 

it follows for the homogeneous model Ai = ... = Am = that 

^ Q^(A; x)^= {-Y^Qu.{-\ -l-x;xf . (3.21) 

Alternatively, we can compute the spectrum of Qoj := O^Qlj^ from the algebraic Bethe 
ansatz similar as before. In terms of the matrix elements of the associated monodromy 
matrices the basic relation we need is 

Qk,iB = - BQk,i + 7^ Qk,i-iA - ^ Qk+i,iD + 3^3' ^ Qk+i,i-iC . (3.22) 

\0k OkOi-iJ di_i dk OkOi-i 



Here the coefficients are the same as in ( p. 10 ). This then leads to the following eigenvalues 
corresponding to Bethe states 

^Q^{X;x)^B{i+)---B{it)\Q) = 

As already previously mentioned for generic inhomogeneity parameters Xm and a suitable 
neighbourhood of a; = (or uj = 00) the Bethe ansatz yields a complete set of eigenstates 



1 23]. This fact now implies the operator equation 



mQUX;x)^ = Q^-^iX + x--x) = Q+_,iX + x)Q-_,iX) . (3.23) 

Therefore, under spin reversal the roles of Q^^, Q^^ are interchanged. These relations match 
the analogous ones derived for the six- vertex model [33, 37, 34 1. 



-10- 



A Q-operator for the twisted XXX model 



4. Fusion hierarchy and complex dimension 

One of the aforementioned advantages of our Q-operator analysis is that the relation be- 
tween Qu;{X;x) and the higher spin transfer matrices t^ is particularly simple allow- 
ing one through analytic continuation to compactly present the information on the entire 
fusion hierarchy. Specializing x — > n G N, it was already pointed out earlier that the 
infinite-dimensional Verma module ( p.8[ ) contains a finite dimensional subrepresentation 
spanned by the vectors {|^)}^Zo which is isomorphic to the SI2 representation tt^""^) 
of spin s = {n — l)/2. The remaining space spanned by {\k)}^^^ can be identified again 
as the Verma module with x = —n. This simple representation theoretic fact trans- 
lates into the following functional relation when splitting the trace over the aforementioned 
subspaces, 

t^^-'H^)=QU^-f,n)-QU^ + f,-n) . (4.1) 

Thus the spectrum of the higher spin transfer matrices takes a particularly simple form 
in terms of the spectrum of Quj{X;x). In contrast the expression from the algebraic Bethe 
ansatz and the fusion relation ( 2.13| ) is more involved. Furthermore, we might analytically 



continue expression ( [4.1[) in the spin variable n setting 

t^(A; x) = Q^{\ -%;x)- Q^{\ + f ; -x) . (4.2) 

The last object combines the information of the entire fusion hierarchy. Notice that in 
(4.1) respectively (4.2) one can safely take the limit to periodic boundary conditions, i.e. 



the following object is well defined 

t(A;x) = lim t^(A;3;) = lim to<^(A - f ;x) - (5^(A + f ; -j;)l . (4.3) 

In this manner one recovers the XXX model with periodic boundary conditions. The 
transfer matrix i(A; x) with "complex dimension" x coincides with the generalized trace 
construction in the context of correlation functions on the infinite lattice. This complex 
dimension occurs as the coefficient of the highest power in the polynomial t{\]x), 

M-l 

t{X; x) = xX^+Y^ tm{x)X^ . (4.4) 

m=0 

In comparison, the analogous result in the context of the six- vertex or XXZ model showed 
the appearance of logarithmic terms; see [p^]. 

4.1 The trace functional: a simple example M = 4, 5*^ = 

It is instructive to verify for a simple example whether the construction (^^) coincides 
with the definition through the trace functional given in |2£]. Setting M = 4 and S"^ = 
we consider a diagonal matrix element of the Q-operator, 

00 

Qc.(A; xZ\::Z = + x-kf{\ + k + If (4.5) 

k=0 
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Here a, = ±1 are the eigenvalues of af acting on the i**^ lattice site with i = 1,2,3,4 
and Qj = 0. Using the formula for the geometric series and analytically continuing the 
result in to afterwards to take the limit w ^ 1 in (^^) we arrive at 

i(A; x)'^l-::Z = ^ + A + A2 + 2x + X (4.6) 

The action of the trace functional Tr^; : U{sl2) <8> C[x] C[x] introduced in |^] (not to be 
mistaken for T^-k^ 7^ Tr^.) on the powers of the Cartan generators is defined through 

^ , sinh(zx) x^ — X 9 7x — lOx^ + Sx^ a , , „, 

^{^ } = ■ \ = ^ + ^ + W^r, ^ + - 4.7 

X smhz 6 360 

Acting now with the trace functional on the monodromy matrix of L-operators we compute 



^ J+^i LX + ^il i A2 + 2Tr{l} A3 + Tr{l} A^ 



t(A;x)2; 



ai...Q;4 
.04 ' 



16 



where the last line is obtained after inserting the values from the expansion ([4.7|). Thus, 
we find agreement with ( [4.31) . To illustrate the generalized transfer matrix of complex 
dimension further we present its eigenvalues in the table below. Specializing x to be an 
integer > one obtains the eigenvalues of each element in the fusion hierarchy. 



p 


t(A;x) 


TT 


2 "2' + ?6 + 'V' A+ A' + 2xA3 + xA^ 


TT 


X x-' x^ , Ax—x'-^ \ 1 Wx—x"^ \2 1 0^ \3 i ™ \4 
6 12 48 + 6 6 A+2XA+XA 





6 + ?6 + ^V' A+ A^ + 2xA3 + xA^ 





X t x-^ i x^ t x-^\ i 2x+x-^ \2 1 0^ \ 3 1 ^ \4 
30 + 12 + 80 + 2 + 2 A+.^XA+XA 


7r/2 


&ix+(i-%i)x''-x'' 1 (4+2i)x-(l+2j)x'=' ^ 1 lOx-x^ \2 | o ^ \3 | ^^4 


7r/2 


8«-(4+8j)x'^+x'' 1 {4-2i)a;'-{l-2i)a;-' | lOx-x^ \^ | X I X X"^ 
48 6 6 



Table 1. Spectrum of the transfer matrix with complex dimension x. 



5. The quantum Wronskian 

The second computational advantage from the Q-operator analysis is of great practical 
importance in the actual computation of the spectra of the Hamiltonian and the transfer 
matrices. Instead of solving the quite intricate Bethe ansatz equations, one can now turn 



the ideology around and rather interpret the relation (O) for x = 1, named the quantum 
Wronskian, as the fundamental identity, 

,^^^^_ u;Q+(A-l)Q-(A)-^-iQ+(A)Q-(A-l) _ ^^^^^ 
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Here we have exploited the factorization ( |3.6D . In terms of the eigenvalues (|3.12| ), (|3.l4 ) 
the above relation incorporates the Bethe ansatz equations above and below the equator 



with respect to the parametrization (2.1 



n J' , I, = n % %^, , n± = M/2 ^ S\ (5.2) 

m=l ?i - Am + i j=l ?j - C,j + i 

and is therefore sufficient to analyze the spectrum. Introducing the elementary symmetric 
polynomials = efc(^^, ^^^) in the Bethe roots 

Q^(A) = EH'e^ A"±-^ (5.3) 

k=0 

the quantum Wronskian ( ^.Ij) becomes the following identity 

/\ X N f ^ \^ '^n-e^M-n-k ~ ^ ^^n-k^M-n-i ,r a\ 

eA/-m Ai,...,Aa/ = > > , n , 5.4 

^-^ ^-^ \m — kl u) — LO ^ 

which is quadratic in the M unknowns e^. Here em(Ai, Am) is the m*^ elementary sym- 
metric polynomial in the inhomogeneity parameters. Furthermore, we use the convention 
= for k < and k > n± = M/2 =p S*^. In contrast the Bethe ansatz equations ( |5.2| ) 
are of order M. The approach based on the quantum Wronskian therefore leads to a 
significant advantage in numerical computations for long spin-chains. 

Note that in the limits u; ^ 0, oo we can easily establish the completeness of the Bethe 
ansatz for generic inhomogeneity parameters by a similar line of argument as it has been 
used in |2^]. Namely, assuming all inhomogeneity parameters to be mutually distinct 
we infer from the quantum Wronskian ( ^.1]) the solutions 



M-n 

n 

j=i " i=i 



a; = oo : Qi(A) = H (A - A^, + 1) and Q" (A) = H (A - A^^.^J 



for any permutation (mi, tum) of the index set {1, 2, M}. Obviously, the number of 
distinct solutions is then (^^) which coincides with the dimension of the associated spin 
sector. For w = the roles of are interchanged. Using the implicit function theorem 
one can then argue that the number of solutions stays the same in the vicinity of the point 
id = oo respectively a; = 0. 

5.1 Special solutions for homogeneous chains of even length and 5^ = 

Let M G 2N and consider the spin sector = 0. Then according to our previous discussion 
and Q~ have the same polynomial degree n = M/2 and in light of ( |3.21| ), ( p.23| ) one 
might expect a simple relationship between them. In fact, based on numerical studies of 
homogeneous spin-chains up to length M = 10 and lo = e**^, (/> G R one confirms that there 
exist 2^/^ solutions to the quantum Wronskian which satisfy 

M G 2N, 5^ = : Q-(A) = (-l)f Q+(-A - 1) . (5.5) 
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Notably, for the mentioned examples M = 2, 4, 6, 8, 10 the eigenvalue of the transfer matrix 
which belongs to the groundstate in the limit a; — > 1 always appears to be among this set 
of special solutions. 

For the numerical investigation it is more convenient to use the second parametrization 
( |2.2| ) of the XXX model, since then the coefficients (not the roots) of the polynomials 

n± 

Qtiu) = i'^^Qti-iu - 1/2) = n - yf) 

are always real numbers. In this parametrization the special relationship ( |5.5| ) becomes 
simply 

MG2N, 5^ = 0: Q-(n) = (-l)f Q+(-n) . (5.6) 

At the moment there is no derivation from first principles for this simplification, however, 
it can be motivated by ( p. 21 ) which states that left and right eigenvectors of the Q-operator 
are related by spin-reversal. As the spin zero sector is invariant under the action of the 
spin-reversal operator it can happen that some left and right eigenvectors of Q coincide 
leading via ( ^.211) to the simplification ( ^.5| ) respectively ( |5.6D . Assuming the latter to 
hold true one can verify it for chains of length M > 10 by inserting this special subset of 
solutions into the Wronskian relation which then simplifies to 

^ ( ^.i^ ^Qtiu - §)Q+(-n - 1) - u:-'Qt{u + §)Q+(-n + §) ^^^^ 



Specializing the spectral parameter to u = Vj' + ^ and u = Vj~ — ^ we now obtain the 
following sets of equations for the Bethe roots Vj~ of this subclass of solutions 

-1 M/2 

(vt + i/^r = — — n - < + + o (5.8) 

to -LO k=l 

and 

M/2 

(vl - = n - vt - i){vl + <), (5.9) 

- w fc=l 

respectively. Since uj lies on the unit circle both equations are equivalent under complex 
conjugation provided the Bethe roots Vj~ are either real or occur in complex conjugate 
pairs. For the mentioned examples this is indeed the case. Dividing these two equations 
yields the familiar Bethe ansatz equations for twisted boundary conditions, 



M 



' — UJ * * 



n 'i I ■ (5-10) 



Thus, we infer that extending the assumption ( |5.6| ) beyond the numerically checked exam- 
ples of spin-chains of length M < 10 is compatible with the Bethe ansatz. The correspond- 
ing eigenvalues of the transfer matrix are of the form 

M/2 

i{u) = i^t{-iu - 1/2) = {uj + Lj-^)u^ + imu^'^"", (5-11) 

m=l 
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i.e. only even powers of the spectral parameter u occur. In addition, the parameters im are 
real and the eigenvalue corresponding to the groundstate in the limit of periodic boundary 
conditions a; — > 1 is distinguished by the fact that all coefficients have the same sign, 
sgnim = sgn(Lij + i^^-*^). We leave a more detailed study of these solutions to future work 
as it involves more extensive numerical calculations. 






M = 10 


TT 

2 


Q+= n5=F0.7769661 ^^-0.3231618 ^^±0.1117312 u2+0.011890969 n =F 0.01189097 
t = ±1.553932 ^^±6.04751 u^±9.74055 u^±7.58483 n2±2. 37584 


vr 
20 


Q+= n5-0.06935158 n^-0.403661 ^3+0.01057672 n2-0.0166721 u - 0.000107935 
t = 1.97538 ?x^°+7.42936 7x^+13.5893 1*^+14.8551 n^+9.33335 ^^+2.59013 


TT 

200 


Q+= 7/5-0.00692881 n^-0.404443 ti3+0.00105731 n2-0.0167203 u - 0.000107938 
t = 1.99975 ti^°+7.49929 n^+13.6758 n^+14.9119 n^+9. 35224 ti2+2. 59241 





Q+= n5-0.404451 ^3-0.0167203 u 

t = 2nio+^ n8+13.6767 u6+14.9125 n'^+9. 35243 n2+2. 59243 



Table 2. Groundstate eigenvalues of the transfer matrix and Q-operator in the spin zero sector for 
various twist parameters. 

5.2 Eigenvalues in the limit of periodic boundary conditions 

Let us make contact with the discussion of Pronko and Stroganov for the XXX model with 



periodic boundary conditions 4> = respectively u; = 1 |24]. Starting from the TQ equation 
on the level of eigenvalues they reported the following quantum Wronskian relation with 
respect to the parametrization (^^), 

Q-{u + i)Q+(n - I) - Q-(u - |)Q+(n + i) = (5.12) 

with the crucial difference that the degree of the second linearly independent solution Q~ 
is now increased by one, 

M-n+l 

2-(-) = ^c?^ n (^-^k)- (5-13) 
k=l 

The degree of the other solution, Q'^, describing the well known Bethe roots above the 
equator remains unchanged, 

n 

Q+(u)= Ui^-vt) ■ (5-14) 

k=l 

We have deliberately denoted their solutions by a different symbol to distinguish them 
from the solutions obtained from our operator construction at quasi-periodic boundary 
conditions. As already pointed out in the introduction the quantum Wronskian ( 5.12| ) has 



a restricted number of solutions which is much smaller than the dimension of the respective 
spin sector fixed by the degree n = M/2 — S^. These solutions must correspond to regular 
Bethe states which yield the highest weight vectors of the various s/2 modules, while the 
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"missing" states are simply descendant states from highest weight vectors which he in a 
different (higher) spin-sector. For instance, in the case of even M the possible number of 
highest weight states in the sector 5^ = is given by (jy^Jg) ~ (m/2-i) 



this number is matched by the solutions to ( 5.12| ); see the table below 



M 


3 


4 


5 


6 


7 


8 


9 


10 




1/2 





1/2 





1/2 





1/2 





dim 


3 


6 


10 


20 


35 


70 


126 


252 


No 


2 


2 


5 


5 


14 


14 


42 


42 



Table 3. Number of solutions to ( 5.12 ) in comparison with the dimension of the spin sector. 



The simplified expression for the transfer matrix in terms of the two linearly indepen- 
dent solutions remains formally the same p3|, however, we remind the reader that the 



degree of Q has changed in comparison with (4.1), 

iiu) = Q~{u + i)Q+{u-i) - Q-{u-i)Q-^{u + i) 



(5.15) 



Naturally, one wonders how the solutions are related to the ones at quasi-periodic 
boundary conditions, Q^{u), when the limit a; ^ 1 is taken. One finds that only a subset 
of the solutions stays finite, the other solutions diverge. In the explicit construction of 
the Q-operator this is due to the fact that the trace over the infinite-dimensional auxiliary 
space does not converge any longer. The number of finite solutions, i.e. those for which 
the limit a; — > 1 is well-defined, approach the solution of Pronko and Stroganov: 



if 



hm |g±(n) 



< oo then 



lini (n) 



(5.16) 



The above relation has been numerically verified for spin-chains up to length M = 10. 
Note that both solutions and Q~ approach in the limit u ^ 1 the same solution 
above the equator. This is to be expected as the degree M — n of Q~ can become smaller 
in the limit of periodic boundary conditions but not greater. At the moment there appears 
to be no Q-operator construction which would yield the other solution Q~ and at the same 
time have the analogous factorization property (3^). The constructions suggested in the 
literature for periodic boundary conditions p5| , p6| , 27] all have degree < M for the spin 
1/2 chain of M sites, while the maximal degree of Q~ is M + 1. 







i{u) 


^,4 _ 3 2 _ J_ 
u 2 " 48 




-i + fn2 + |n4 + 2n6 


,,4 1 4=Fv^„,2 7=f2v^ 
U -1- 2 " 16 




31±Wl3^7±8v'l3^2+9^4 + 2n6 


" ^ ^ ^ 2v^ 16 


,,3 _|_ JL 4- ^ 
^ ^ 12 =^ 4v^ 


=F VSu + f + + 2u^ 



Table 4. Solutions to the quantum Wronskian ( 5.12 ) for M = 6, S^= and the corresponding 
eigenvalues of the transfer matrix. 
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6. Conclusions 



In this work we have presented the isotropic hmit of a previous Q-operator construction 
for the XXZ model 33, 34 1 in order to discuss the XXX model with quasi-periodic 
boundary conditions. The motivation for this discussion has been twofold. On the one 
hand this construction enables one to formulate an analytic continuation of the fusion 
hierarchy to complex dimension as it has been recently used in the description of correlation 



functions in form of a trace functional |2£]. In this context it should be noted that previous 



constructions of Q-operators for the XXX model |25, 27 1 have always been for periodic 
boundary conditions where an analogous formulation does not exist. This is due to the 
fact that the trace over an infinite-dimensional auxiliary space has to be taken whose 
convergence is not necessarily guaranteed. Moreover, due to the sl2 symmetry the set of 
solutions to the Bethe ansatz equations is reduced (i.e. only the highest weight states in 
each s/2 module are proper Bethe states), whence certain functional relations such as the 
quantum Wronskian for periodic boundary conditions do not yield the complete set of 
eigenvalues; compare with table 3. 

This provided additional motivation for investigating a Q-operator for the twisted 
XXX model. Via this construction one is lead to a quantum Wronskian for quasi-periodic 



boundary conditions (see (5.1) in the text), which now yields the complete set of Bethe 
states and eigenvalues of the transfer matrix. Our derivation relied on previous algebraic 
Bethe ansatz results for the Q-operator of the XXZ model |^^. As emphasized in the text 
the quantum Wronskian has a simpler structure than the Bethe ansatz equations and based 
on numerical computations we found special solutions for spin-chains of even length and 
vanishing total spin satisfying more fundamental identities. For instance the Bethe roots 
of the aforementioned subset of solutions obey the set of equations. 



M 



LO 



M/2 

k=i 



and are either real or occur in complex conjugate pairs; see the discussion in section 5.2. 
Among these special solutions is the eigenvalue which corresponds to the groundstate in the 
limit of periodic boundary conditions and has real Bethe roots. The present numerical data 
only include chains up to length M = 10 and further investigation is needed to see whether 
they persist for longer chains. This is particular important in order to make contact with 
the thermodynamic Bethe ansatz and the string hypothesis [10| |38| [16|. As it has been 
discussed in the literature there might be a critical length beyond which certain solutions 



cease to exist, see e.g. |10] |39]. We leave this problem of a more extensive numerical study 
to future work. 
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